Collective orbital excitations in orbitally ordered YVO3 and H0VO3 
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We study orbital excitations in the optical absorption spectra of YVO3 and H0VO3. We focus on 
an orbital absorption band observed at 0.4 eV for polarization E \\ c. This feature is only observed in 
the intermediate, monoclinic phase. By comparison with the local crystal-field excitations in VOC1 
and with recent theoretical predictions for the crystal-field levels we show that this absorption band 
cannot be interpreted in terms of a local crystal-field excitation. We discuss a microscopic model 
which attributes this absorption band to the exchange of two orbitals on adjacent sites, i.e., to the 
direct excitation of two orbitons. This model is strongly supported by the observed dependence on 
polarization and temperature. Moreover, the calculated spectral weight is in good agreement with 
the experimental result. 



I. INTRODUCTION 

In strongly correlated transition-metal oxides, orbital 
interactions play a key role in many intriguing phenom- 
ena such as the colossal mag netores istance or the effective 
reduction of dimensionality.'^ 3 ^ Orbitals on different 
sites interact with each other^l via the collective Jahn- 
Teller effect, i.e., the coupling to the lattice, and via ex- 
change interactions, which are governed by the antisym- 
metrization of the total wave function including both the 
orbital and the spin part. These interactions can result 
in coupled long-range spin and orbital order. If the cou- 
pling to the lattice is dominant, the excitations are well 
described by "local" crystal-field (CF) excitations,* 7 * 8 * 9 ' 
where "local" means that the excitation can be treated 
as a change of the orbital occupation on a single site, i.e., 
the dispersion is negligible. In the opposite case of dom- 
inant exchange interactions, one expects novel collective 
elementary excitations, namely orbital waves (orbitons) 
with a significant dispersion ^reflecting the propagation 
of the excited state. Thus orbitons are analogous to spin 
waves - propagating spin flips - in a magnetically or- 
dered state. Orbitons are expected to reveal the funda- 
mental orbital interactions responsible for the interesting 
physical properties. In the quest for the experimental 
observation of orbitons, the central experimental task is 
to demonstrate that the orbital exchange interactions are 
essential for the elementary excitations. If this is the case, 
the excitations cannot be described in terms of single-site 
physics, and we will use the term "orbiton" . 

The first claim for the observation of orbitons was 
based on Raman data of LaMn03,12 but the relevant 
features later have been explained in terms of multi- 
phononsP^In fact, in the manganites the orbital degree 
of freedom is connected with e g electrons, for which the 
coupling to the lattice is strong in an octahedral envi- 



ronment. The vanadates i?V03 with two electrons occu- 
pying t2 g orbitals may be considered as more promising 
candidates ! 13 ' 14 ' 1 ^ Recently, the observation of orbitons 
in Raman data of -RVO3 (R=Y, La, Nd) ha s been claimed 
at 43 and 62meV by Miy asaka et a?! 17 ' 18 ' and at 45 and 
84 meV by Sugai et al. ™ but the proposed orbitons are 
hard to discriminate from (multi-)phonons and magnons, 
and the assignment is controversial .liZHSllSJ Thus, an ex- 
perimental proof for the existence of orbitons is still lack- 
ing. 

Here, we report on the observation of orbital excita- 
tions in the optical conductivity a(u>) of orbitally ordered 
YVO3 and H0VO3. We focus on an absorption feature 
observed for E || c at about 0.4 eV, well above the range of 




(a) T < 77K (Pbnm) (b) 77K < T < 200K (P2/b) 



FIG. 1: Orbital and spin ordering patterns in YVO3 
for ( a) T s < 77 K and (b) T s < T < Too = 
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pied by one electron on each site (not shown) . The occupation 
of xz and yz orbitals is C-type below Ts and G-type above 
(see main text for more details). Spin order (SO) is lost at 
Tjv = 116 K. The coordinates of the Pbnm and P2\/b crystal 
systems are given in the upper left corner of each figure. 
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phonons and magnons and well below the Mott- Hubbard 
gap. A comparison with the local CF excitations i n the 
3d 2 system VOC1 and with recent theoretical result 
for the CF levels shows that this feature is hard to recon- 
cile with a local CF scenario, in particular as far as the 
energy, the polarization and temperature dependence are 
concerned. We discuss the microscopic exchange process 
and conclude that the feature at 0.4 eV reflects the ex- 
change of two orbitals on neighboring sites, i.e., the direct 
excitation of two orbitons. 

The vanadates -RVO3 with R = Y and Ho exhibit 
an orthorhombic crystal structure (Pbnm) at room 
temperatureP^^' The undoped compounds represent 
Mott-Hubbard insulators with two localized electrons 
in the 3d shell of each V 3+ ion. A crystal field of 
predominantly octahedral symmetry yields a splitting 
of the 3c? states into a lower-lying, triply degenerate 
tig level and a doubly degenerate e g level. A de- 
tailed analysis of the structure reveals that the degen- 
eracy of these levels is fully lifted by an orthorho mbic 
distortion of the VOe octahedra (Dih symmetry)) 20 * 21 ! 
giving rise to a splitting of the t 2g manifold into xy, 
xz, and yz orbitals. In YVO3, a low-temperature 
orthorhombic phase (Pbnm) with G-type spin order 
(SO) and C-type orbital order (00) w as found below 

Ta =77K pll6l20l21l22l23l24l25l26l27b8l29l30l ^ the xy QJ ._ 

bital is occupied at each V site, whereas the occupation of 
xz and yz orbitals alternates within the ab plane (see Fig. 
[IJ . At T5 = 77 K a first-order structural phase transition 
to an intermediate, monoclinic phase with P2i/b symme- 
try has been observed. This monoclinic phase shows two 
different vanadium sites, V(l) and V(2), which alternate 
along the c axis. Therefore, the mirror symmetry per- 
pendicular to the ab plane is broken in the intermediate 
phase. The spin ordering pattern changes from G-type 
below Ts to C-type above. Long-range magnetic order 
is lost at T/v = 116K. The structural phase transition 
from the monoclinic phase to the orthorhombic room- 
temperature phase is observed at Too ~ 200 K, which 
has been interpreted as the long-range orbital ordering 
temperature™ However, synchrotron x-ray diffraction 
data give evidence for the presence of OO up to about 
300K.^SlFor the intermediate monoclinic phase, it has 
been discussed controversially whether the physics has to 
be described in terms of 'classical' orbital order or quan- 
tum orbital fluctuations: It has been proposed that the 
intermediate phase of YVO3 represents the first realiza- 
tion of a one-dimensional orbital liquid and of an orbital 
Peierls phase with V(l)-V(2) orbital dimers^mMMM 
This h as bee n challenged by LDA+U and LDA+DMFT 
studiesp^H which for YVO3 find orbital order and that 
at least below 300 K orbital quantum fluctuations are 
suppressed in YV0 3 by a sizeable ligand-field splitting. 
The orbital ordering pattern in the monoclinic phase has 
been reported as G-type based on, e.g., resonant x-ray 
diffraction,^ an analysis of the V-0 bond lengths,^ or 
LDA+U calculations.^ In comparison to C-type OO, the 
orbitals of every second layer along c are shifted along x, 
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FIG. 2: (color online) Optical conductivity of YVO3 in the 
low-temperature orthorhombic phase at T— 10 K for E \\ a, b, 
and c (top panel) and -ln(T(o>)) of VOC1 at 20 K for E \\ a 
and b (in Pmmn; bottom panel). Three different orbital exci- 
tations within the £29 subshell are sketched in the lower panel. 

thus xz and yz alternate along x, y, and z for G-tye OO 
(see Fig. [I]). According to Hartree-Fock calculations, the 
size of GdFe03-type distortions is decisive for the choice 
between G-type and C-type orbital order!^ A recent 
LDA+DMFT study by De Raychaudhury et al. finds 
that the OO pattern is intermediate between C-type and 
G-type due to the GdFe03 distortions, almost C-type 
in the intermediate phased However, Noguchi et alW* 
claim that their synchrotron x-ray diffraction data are 
fitted best by G-type OO, but that the partial occupation 
of other orbitals is possible. The compound H0VO3 be- 
haves very similar to YVO3 with slightly different phase- 
transition temperatures of Too ~ 188 K, IV = 114K, and 
Ts ~ 40 K.I2H2SI Detailed neutron and synchrotron x-ray 
scattering experiments were unable to establish the pres- 
ence of orbital fluctuations in H0VO3P2I 

II. EXPERIMENT 

Single crystals of .RVO3 with R = Y and Ho have been 
grown by the traveling-solvent floating-zone method.^ 
The purity, stoichiometry and single-phase structure of 
the crystals was checked by x-ray diffraction and ther- 
mogravimetry. Typical crystal dimensions are a few mm 
along all three crystallographic axes. The optical con- 
ductivity o(uj) of YVO3 was determined by measuring 
both the transmittance and the reflectance^ between 
0.06 and 1.9 eV using a Fourier spectrometer. The mea- 
surements have been performed using linearly polarized 
light with the electric field parallel to the orthorhombic 
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axes, i.e., E \\ a, 6, and c (see top panel of Fig. [2|. For 
convenience, we use the same set of axes at all temper- 
atures, i.e., we neglect the monoclinic distortion of the 
structure. This is justified because the monoclinic an- 
gle a = 89.98° is very close to 90° The reflectance was 
measured on samples with a thickness of d > 2 mm in or- 
der to avoid backside reflections. The transmittance data 
were collected on a series of crystals with different thick- 
ness (100 fjxn < d < 500 ^m), which were polished on both 
sides. 

Single crystals of VOC1 have been grown by the 
chemical-vapor transport technique. The purity of the 
crystals was checked by x-ray powder diffraction. Typi- 
cal crystal dimensions are a few mm 2 in the ab plane and 
10-100 p,m along the c axis. In the case of VOC1 and 
H0VO3, we have measured the transmittance only. The 
transmittance T(w) is a very sensitive probe for the de- 
termination of weakly infrared-active excitations below 
the gap of these Mott-Hubbard insulators, where the re- 
flectance is nearly constant and featureless. Therefore, 
the absorption coefficient a(uS) cx — ln(T(w))/d can be 
used equivalcntly to u(uj) for the determination of weak 
orbital excitations.^ 



III. RESULTS AND DISCUSSION 

The top panel of Fig. [2] shows the optical conductiv- 
ity cr(u>) of YVO3 at 10 K in the transparent window 
of the Mott-Hubbard insulator, i.e., above the phonon 
range and below the electronic interband excitations. 
The lowest electronic transition corresponds to an ex- 
citation across the Mott-Hubbard gap, i.e., the trans- 
fer of one electron from a 3d 2 V 3+ site to another one, 
3d 2 3d 2 — ► 3d 1 3d 3 . In the optical conductivity, the low- 
est electronic excitation is observed at 1.8 eV |25l38j . 
and a value of 1.6 eV has been reported for the Mott- 
Hubbard gap based on a combination of ellipsometry 
and LSDA+I7 calculations!^ This agrees with the steep 
increase of <j{uj) in our data (see top panel of Fig. |2j. 
However, the transmittance reveals that the very onset 
of excitations across the Mott-Hubbard gap is somewhat 
lower. The precise onset is obscured by the superposi- 
tion of spin-forbidden orbital excitations between about 
1.0 and 1.5eV (see below). Charge-transfer excitations, 
involving the transfer of an electron between V and O 
ions (d 2 -^d 3 L), are located above about 4 eV.^ Absorp- 
tion features below the Mott-Hubbard gap have to be at- 
tributed to phonons, magnons, excitons, orbital excita- 
tions, or to localized carriers trapped by impurities. We 
exclude the latter for a number of reasons: (i) the spectra 
of different samples of YVO3 are identical, (ii) the spectra 
of YVO3 and H0VO3 are very similar (see below), (iii) the 
polarization and temperature dependence (see below), 
and (iv) the DC resistivity, which is very large, 3T0 6 ficm 
at 200 K. In earlier experiments, an activation energy of 
A act = 0.25eV was obtained for polycrystalline samples 
from resistivity data for 180 K< T < 300 K, 39 correspond- 
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FIG. 3: (color online) (a) Resistivity p of YVO3 from 
T = 200K up to 570 K. (b) In (p) vs. 1/T (left axis) and 
ln(p/T5) V s. 1/T (right axis). Dashed lines show fits for 
activated behavior (left axis) and within the non-adiabatic 
small-polaron model (right). 



ing to an optical gap of 2A act =0.5eV. In Fig. [3^a) we 
plot p(T) of our single crystals between 200 and 570 K. 
Note that p(T) is about an order of magnitude larger 
than reported in Ref. |39l The Arrhenius plot in Fig. (3^b) 
suggests 2A act ps 0.8 eV, but the data are not very well de- 
scribed by simple activated behavior. The data are best 
described by the non-adiabatic small-polaron modelj^ 
predicting p= CT§ exp (A act /(fc B T)) (see Fig. gb)). 
Within this model we obtain 2A act = 0.78 eV, in agree- 
ment with both the estimate based on simple activated 
behavior and in particular with the onset of absorption 
observed in a(uj) (see top panel of Fig. |2j). This value is 
significantly larger than reported in Ref. 1391 demonstrat- 
ing the high quality of our samples. Remarkably, DFT- 
PIRG calculations 3 ^ predict an indirect gap at 0.7 eV, in 
excellent agreement with our data. Altogether, we have 
strong evidence that all absorption features observed in 
a{uj) below 0.8 eV arise from phonons, magnons, exci- 
tons, or orbital excitations. 

The optical conductivity of YVO3 below 0.9 eV is given 
in Fig. [4] The absorption due to phonons gives rise 
to the steep increase of <r(oS) below about 80meV, in 
agreement with the data reported in Ref. |25l Thus, 
weak two-phonon features can be observed up to about 
160 meV, whereas three-phonon absorption is expected to 
be still much weaker. A contribution to cr(u>) from spin 
waves may arise in the form of two-magnon-plus-phonon 
absorption. 3D However, in YVO 3 the spin- wave disper- 
sion does not exceed 40meV!221 Thus, a possible two- 
magnon-plus-phonon contribution is expected to peak 
below 0.2 eV and clearly is not related to the features 
observed above 0.2 eV. 

As far as excitons are concerned, we have to distin- 
guish between strongly and weakly bound excitons. In 
a Mott-Hubbard insulator, an exciton is a bound state 
of an electron in the upper Hubbard band (i.e., a dou- 
ble occupancy) and a hole in the lower Hubbard band 
(an "empty site"). In the case of a weakly bound exci- 
ton, the electron and the hole occupy distinct sites, e.g., 
nearest-neighbor sites. The binding may result from the 
nearest-neighbor Coulomb attraction. Recently, an exci- 
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tonic resonance was reported^ in the orbitally ordered 
3d 1 Mott-Hubbard insulator YTi0 3 at 1.95 eV. We con- 
sider it as unrealistic to assume a binding energy large 
enough to pull such a weakly bound exciton far below 
leV in YVO3. Moreover, the spectral weight of the ex- 
citon in YTiC>3 is more than two orders of magnitude 
larger than the spectral weight between 0.1 and 0.7 eV in 
YVO3. 



Orbital excitations 

The physics is different for a strongly bound exciton, 
in which case the electron and the hole share the 3c? shell 
of the same site. This actually corresponds to an or- 
bital or d-d excitation (or to a spin excitation, which has 
been ruled out above). These excitations may occur far 
below the gap because the binding energy and the gap 
are both of the order of the on-site Coulomb repulsion 
U. In other words, double occupancy is omitted if the 
electron and the hole share the same site, thus one does 
not have to pay the energy U. Orbital excitations in 
the form of local crystal-field excitations are a common 
feature in many transition-metal compounds in the con- 
sidered frequency range EEEl The central question to us 
is whether the observed orbital excitations are such lo- 
cal crystal-field excitations or propagating orbitons, re- 
flecting the importance of orbital exchange interactions. 
The orbiton dispersion for orbitally-ordered vanadates 
has been inves tigated theoretically.^ Predictions for Ra- 
man scatterin g* 13 ! 17 ! 18 ! anc j inelastic neutron scattering^ 
have been discussed, but contributions to <r(u>) for the 
case of dominant exchange interactions have not been 
considered thus far. Moreover, a quantitative description 
requires that both the exchange interactions and the cou- 
pling to the lattice are treated on the same footing,S2EH 
but up to now such calculations have not been reported 
for the vanadates. Therefore, we first compare our results 
with the expectations for local crystal-field excitations of 
3d 2 V 3+ ions. We use VOC1 as a typical example for the 
absorption spectrum of V 3+ ions in a predominantly oc- 
tahedral, but distorted environment (see bottom panel of 
Fig. [2j) . In the sister compound TiOCl, the orbital excita- 
tions are v ery w ell described in terms of local crystal-field 
excitations P3S 

YVO3 shows inversion symmetry on the V sites, thus 
local CF excitations are not infrared-active due to the 
parity selection rule. However, they become weakly al- 
lowed by the simultaneous excitation of a symmetry- 
breaking phonon! 7 ^ 8 ^^ VO C1 (P mmn) there is no inver- 
sion symmetry on the V sit d 46 l 47 l and the CF excitations 
are weakly allowed without the additional excitation of 
a phonon. 

For the sake of simplicity, we start from an undistorted 
octahedral crystal field, the effect of a lower symmetry 
will be discussed below. The ground state of a d 2 system 
in an octahedral field is the nine-fold degenerate 3 Ti level 
with total spin S = 1, in which two electrons with paral- 
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FIG. 4: (color online) Temperature dependence of o(lo) of 
YVO3 in the mid-infrared range for E \\ a (top panel), E \\ b 
(middle), and E \\ c (bottom). The features at 0.4 eV for E \\ c 
and 0.55 eV for E \\ a constitute our main experimental result. 

lei spins occupy the ti g level, t^ g The splitting between 
t2 g and e„ levels typically amounts to > 2eV for O 2- 
ligands. 7 8 For VOC1, this splitting is reduced due to the 
smaller ligand strength of the CI 1- ions and can be identi- 
fied with the feature observed around 1.7 eV (see bottom 
panel of Fig. |2j. For YVO3, it is reasonable to assume 
that the ti g -e g splitting is larger than the Mott-Hubbard 
gap. In the following, we focus on the excitations within 
the t2 9 shell, which are located at lower energies. 

In an octahedral field, the spin-flip excitation from the 
3 Ti ground state with S — 1 to the five-fold degenerate 
5 = state { X T 2 , 1 E; t\ g ) occurs at 2J H . The Hund ex- 
change J H ss 0.7 eV (Ref.[25) is hardly screened in a solid, 
therefore this excitation is obs erved at very similar ener- 
gies in different V compounds] 50 * 51 !^ irrespective of the 
crystal structure. Typical examples are the weak, sharp 
features observed between 1.1 and 1.3 eV in VOC1 (see 
bottom panel of Fig. |2j. These spin-flip excitations are 
very weak due to the spin-selection rule, they become 
weakly allowed by spin-orbit coupling or by the simulta- 
neous excitation of a magnonPl Such spin-flip bands often 
are very sharp because the orbital occupation stays the 
same, thus the coupling to the lattice is only weak (the 
large width of spin-allowed excitations is attributed to vi- 
bronic Franck-Condon sidebands). Also in YVO3, sharp 
features are observed between 1.1 and 1.4 eV which can 
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be attributed to spin-forbidden excitations. The oscil- 
lator strength is clearly enhanced compared with VOC1. 
Presumably, this is due to the overlap with the onset of 
excitations across the Mott-Hubbard gap. Mixing the 
two kinds of excitations will transfer some weight to the 
orbital bands, a process called "intensity stealing". 

In a local crystal-field scenario, all features observed 
significantly below 1 eV have to be interpreted as spin- 
conserving transitions within the shell, i.e, both the 
ground state and the excited state show S =1. Both in 
YVO3 and in VOC1, the symmetry on the V sites is lower 
than tetragonal, thus the t2 g level is split into three dis- 
tinct orbitals. For two electrons with parallel spins, there 
are three distinct energy levels, each showing three-fold 
spin degeneracy In strongly distorted VOC1 we observe 
the corresponding excitations at 0.3-0.4eV (see Fig. |2j. 
For YVO3, one expects lower excitation energies because 
the distortions away from an octahedral environment are 
smaller than in VOC1. Indeed, recent first-principles^ 
and LDA+DMFT calculations^ predict intra-i2g excita- 
tions in YVO3 in the range of 0.10 - 0.20 eV and 0.06 - 
0.24 eV, respectively. Our data of cr(w) of YVO3 show 
an absorption band centered around 0.20 - 0.25 eV at all 
temperatures and for all polarization directions (see Fig. 
[3]), in agreement with these expectations. 

Our main experimental result is the observation of the 
two remaining features which cannot be explained within 
a local crystal-field scenario, namely the peaks at about 
0.4eV for E \\ c and at 0.55 eV for E \\ a (see Figs. [I] and 
[5]). As follows from the discussion above, these features 
are hard to reconcile with a local crystal-field scenario: 
the energy is too low for a spin-flip transition and too 
high for a spin-conserving intra-i2g transition. Moreover, 
the strong polarization dependence of both peaks is en- 
tirely unexpected for a phonon-assisted crystal-field ex- 
citation. Since the phonon polarization is arbitrary, one 
does not expect strict polarization selection rules. 9 Even 
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FIG. 5: (color online) Spectral weight and line shape of the 
peak at 0.4 eV for E\\c in YVO3. According to our micro- 
scopic model, the two-orbiton excitation is forbidden below 
Ts = 77 K. Therefore, the data for 4K have been subtracted 
for each temperature as an estimate for the background of, 
e.g., multi-phonon absorption. The inset depicts the plasma 
frequency lu p (see Eq. |8| for uj\ = 0.20 eV and w hi = 0.68 eV. 
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FIG. 6: (color online) Top: Sketch of the local energy levels 
and of the hopping processes between neighboring V sites for 
the two-orbiton excitation in the monoclinic phase of KVO3 
with an electric field E applied parallel to the c axis. Only 
xz and yz orbitals are considered for the exchange because 
the lower-lying xy orbital is always occupied by one electron. 
The spin is omitted since we consider fully polarized elec- 
trons. See text for more details. Bottom: For two electrons 
with parallel spins, there are six possibilities to occupy the 
xz and yz orbitals on sites V(l) and V(2): (a) orbital ground 
state; (b) two-orbiton excited state; (c) & (d) one-orbiton ex- 
cited states; (i), (ii) states in the upper Hubbard band with 
a doubly occupied V site. 



in VOC1, where CF excitations are directly allowed, the 
absorption is very similar for the two polarization direc- 
tions. Finally, also the pronounced temperature depen- 
dence is unexpected for a local CF excitation. For E\\c, 
the spectral weight around 0.4 eV is independent of tem- 
perature for T > Tqo = 200 K, increases with decreas- 
ing temperature below Too, an d abruptly disappears at 
T S = 77K (see inset of Fig. [5]). For a local CF excita- 
tion, both the spectral weight and the energy in prin- 
ciple may change across a structural phase transition. 
However, LDA+DMFT calculations^ predict that the 
intra-^2g excitation energies change by less than 40 meV 
across Ts, thus the abrupt change of cr(u>) at Ts cannot 
be explained. 



Orbitons 

Thus far we have neglected the exchange interactions 
between orbitals on neighboring sites, which change the 
character of the excitations from "local" CF excitations 
to propagating orbitons. Here, we consider two differ- 
ent processes for the excitation of orbitons.^ First, the 
excitation process itself may be based on the exchange 
of orbitals between adjacent sites (see below). Second, 
an orbiton may be excited locally by flipping an orbital 
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on a single site, e.g., from xz to yz, with subsequent 
propagation of the orbital flip. In the latter case, the ex- 
citation process is as discussed above in the local crystal- 
field limit, i.e., it requires the simultaneous excitation 
of a symmetry-breaking phonon to obtain a finite dipolc 
moment for the local orbital flip. The difference arises 
from the propagation of the orbital flip, which is sim- 
ilar to a spin flip which translates into a spin wave or 
magnon in a system with long-range magnetic order. For 
such a propagating orbital flip we have to take into ac- 
count the dispersion. Due to momentum conservation, 
a(uj) is sensitive to excitations with Ak to t = 0, where 
Ak to t = kphonon+korbiton- Since the symmetry-breaking 
phonon may have any momentum, o~(u>) reflects orbiton 
contributions from the entire Brillouin zone. However, 
without more detailed theoretical predictions about the 
line shape or width, this one-orbiton-plus-phonon peak 
cannot be distinguished from the broad vibronic Franck- 
Condon peak expected in the local CF limit. Thus, we 
can attribute the absorption feature around 0.2 eV to a 
single orbital flip, but we cannot decide whether this flip 
is only local or propagates due to orbital exchange inter- 
actions. 

In the other case mentioned above, the excitation pro- 
cess is based on the exchange of orbitals between adja- 
cent sites. Here, we primarily focus on the peak for E\\c 
at 0.4 eV. Our interpretation of this feature as a two- 
orbiton excitation naturally explains its energy, the po- 
larization dependence, and the pronounced temperature 
dependence, i.e., it resolves the three issues mentioned 
above. For the exchange along the c axis, we consider 
only the dominant processes and neglect the rotation and 
tilt distortions of the octahedra. In this case, hopping 
preserves the type of orbital, i.e., the only finite hopping 
processes are xz{a) <-> xz(/3) and yz{a) <-» yz(/3), where 
a and [3 denote neighboring V(l) and V(2) sites, respec- 
tively (see Fig. [6]). Note that hopping in z direction is 
zero for the xy orbital, which is the lowest orbital on all 
sites and thus occupied by one electron. Following the 
analysis of synchrotron x-ray data by Noguchi et aZ.p^we 
consider G-type orbital order for the intermediate phase 
(see Fig. i.e., the second electron per site occupies xz 
on V(l) and yz on V(2) in the ground state. For this 
second electron per site we consider the fermionic Hamil- 
tonian Hp 

Hp = i?Fo + i(4i c /3o + c lo c /3i + h.c.) (la) 
H F0 = ((e T + A T )4 1 c rl +£ T 4 c T0 

t£{q,(3) 

+ c/'4 lCTl 4 c r0 ) (ib) 

where t denotes the hopping matrix element, c ri (c ri ) 
creates (annihilates) an electron in orbital i = or 1 on 
site r £ {a,f}}, U' = U — 3Jh denotes the Coulomb re- 
pulsion for two electrons on the same site but in different 
orbitals, and the energies e T and A T are illustrated in 
Fig. U No spin appears since we consider parallel spins 



due to Hund's coupling to the electron in the xy orbital. 
In the orbital ground state, the lower levels (denoted 0) 
are occupied by one electron, the upper levels (denoted 
1) are empty. The creation of an orbiton at site r cor- 
responds to the excitation of an electron from to 1 at 
r, which requires the energy A T . We introduce bosonic 
orbiton creation operators b\ with 

4 == 4ic.ro- (2) 

The annihilation operators are the hermitean conjugate 
ones. Note that the orbitons are hardcore bosons since 
at maximum there can be only one at each site. 

In second order in t the orbital ground state of Hp 
(state (a) in Fig. [6]) is linked via the intermediate states 
(i) or (ii) to state (b) with two excited orbitons, one at 
site a and one at site /3. Note that states (c) and (d) 
corresponding to the excitation of a single orbiton can- 
not be reached from the ground state in the considered 
symmetry.^1 In contrast, the two-orbiton excitation can 
account for the observed energy, polarization and tem- 
perature dependence. For the energy of a two-orbiton 
excitation, one roughly expects twice the energy of a 
one-orbiton excitation, neglecting orbiton-orbiton inter- 
actions and kinetic effects. As discussed above, single 
orbital excitations are observed around 0.2 eV, hence the 
energy of the feature at 0.4 eV is well described by a two- 
orbiton interpretation. Due to Hund's coupling with the 
electron in the low-lying xy orbital, this two-orbiton ex- 
citation requires parallel alignment of the spins on the 
considered bond, which in YVO3 is only present in the 
intermediate phase and only along the c axis, in agree- 
ment with the observed polarization and temperature de- 
pendence. Moreover, the two-orbiton excitation requires 
hopping from, e.g., xz on V(l) to xz on V(2), thus the 
spectral weight becomes zero if both orbitals are occu- 
pied in the ground state. This is the case for the C-type 
OO observed below Tg, explaining the abrupt drop of 
the intensity at Tg. For the intermediate phase, we have 
assumed pure G-type OO for simplicity, neglecting ad- 
mixtures of C-type (see discussion in the introduction). 
Note that deviations from G-type OO as claimed recently 
on the basis of LDA+DMFT calculations^ will only af- 
fect the precise value of the spectral weight. In order to 
understand the temperature dependence at higher tem- 
peratures, we now address the dipole-selection rule and 
calculate the spectral weight of this two-orbiton excita- 
tion. 

In analogy to the derivation of Heisenberg exchange, 
we derive an effective Hamiltonian iJ or b in terms of or- 
biton creation and annihilation operators for the mixture 
of states (a) and (b) 

H rb = J(blbl + h.c.) + const. (3) 

In this effective description no virtual double occupancies 
appear. To do so we assume that the on-site Coulomb 
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interaction is the largest energy 

U' > \e a + A a -ef3\ (4a) 
U' > \ep+A p -e a \. (4b) 

This assumption is certainly met in YVO3 where U' = 
U — 3Jh ~ 2 eV and all other energies are of the order of 
fractions of an eV. We find (see Appendix) 

2t 2 U' 

(U') 2 - (6 - ecEf ' 1 ' 

where e is the elementary charge, c the distance between 
the two sites, E denotes an electric field applied parallel 
to the bond along the c axis, and 

^ = 4-4 + A ^ a = (4+4-4-4)/2. (6) 

The two-orbiton excitation is dipole allowed if 
8H mh /dE^0. We find 

At 2 

dJ/dE^-^Sec. (7) 

In the denominator we neglected 6 because it is much 
smaller than U'. 

In the presence of a mirror plane between the two V 
sites, i.e., e° = e° a and = A a , the two-orbiton ex- 
citation does not carry a dipole moment and does not 
contribute to a(u>). However, such a mirror plane is 
present only below Ts = 77 K and above Too = 200 K, 
but the symmetry is broken in the intermediate phase 
with two distinct V sites. The situation is similar to the 
case of two-magnon absorption discussed by Lorenzana 
and SawatzkyPH They demonstrated that two-magnon 
absorption becomes weakly infrared active if the mirror 
symmetry on the bond is broken by the simultaneous ex- 
citation of a phonon. In the present case, the symmetry 
is already broken without a phonon. We conclude that 
the excitation of two orbitons is directly infrared active 
for E || c in the intermediate, monoclinic phase of YVO3. 

The absolute value of dJ/dE depends on the crystal- 
field levels via 5. From the first-principles study by 
Solovyev^and the LDA+DMFT calculation by De Ray- 
chaudhury et oi.™ we obtain 6 — 20 and 33.5 meV, re- 
spectively. However, each individual level e\ certainly is 
known only to within 50 meV, thus the error of 6 is about 
100 meV. For U' = 2 eV, S = 20 - 100 meV, and t = 100 - 
150 meV (Ref.ESJ) we obtain OJ/dE = 0.1 - l.M0- 3 ec. 
The spectral weight usually is expressed in terms of the 
plasma frequency u> p , 

ai{u))du = — uj 2 , (8) 

where the frequency range of interest is defined by lu\ 
and Whi) an d £0 denotes the dielectric constant of vac- 
uum. The spectral weight of the two-orbiton excitation 
at hx>2o = 0.4 eV is given by 

^ = ^\9J/dE\\ (9) 



where V = 56 A 3 is the volume per site. Finally we obtain 
hujp = 0.6 - 6.8 me V. This has to be compared with the 
experimental result for the spectral weight, for which we 
choose hcu\ Q — 0.20 eV and Whi = 0.68 eV. In order to sepa- 
rate the two-orbiton contribution from a background of, 
e.g., multi-phonon absorption, the integration for each 
temperature T is performed over cr 1 (T, w) — tri(4K,a;), 
since the two-orbiton excitation is not dipole allowed at 
4 K (see Fig. [5h . At 80 K this yields hu c ^ = 16 meV, 
about a factor of 2 - 26 larger than the calculated re- 
sult. We emphasize that the quantitative prediction of 
the spectral weight is a challenging task. Note that the 
theoretical estimate of oj p of the phonon-assisted two- 
magnon absorption in the cuprates on the basis of a sim- 
il ar pe rturbation expansion was off by a factor of 4 - 
7 mm Therefore we consider this result as a clear sup- 
port for our interpretation. 

The temperature dependence of u> p is plotted in the 
inset of Fig. [5] Upon warming above Too = 200 K, the 
mirror plane is restored and the direct contribution to 
<j(uj) is suppressed. However, the spectral weight above 
Too is larger than at 10 K. This may either arise from a 
weak, phonon-assisted two-orbiton contribution or from 
thermal broadening of the electronic gap and of the ab- 
sorption band around 0.2 eV (see Fig. |5|. In the case 
of phonon-assisted two-magnon absorption proposed by 
Lorenzana and Sawatzky, 41 <j(oj) in low-dimensional an- 
tiferromagnets is hardly affected at the magnetic order- 
ing temperature Tv, because the spin-spin correlation 
length remains large above Tv. In YVO3, measurements 
of the specific healP^l show that only of the order of 
10% of the expected entropy are released in the phase 
transition at Too, which indicates strong fluctuations. 
Therefore, a finite contribution of phonon-assisted two- 
orbiton excitations is possible above Too- I n contrast, 
the spectral weight abruptly vanishes upon cooling be- 
low Ts = 77 K, because the mirror symmetry is restored 
and because both the orbital and the magnetic ordering 
patterns change (see Fig. [IJ. The reduction of Hlj p from 
16 to lOmeV observed between Ts and Too can tenta- 
tively be attributed both to a reduction of ferromagnetic 
correlations between nearest neighbors above T/v and to 
rather small changes of the orbital occupation or of the 
crystal-field levels. For 5 = 20 - 100 meV, a reduction of 
Up by a factor of 1.6 corresponds to changes of the four 
individual CF levels of only about 4-20 meV. 

The line shape in principle may serve as a key feature 
to test our interpretation. However, this requires to take 
into account the orbital exchange interactions and the 
coupling to the lattice on the same footingPSl On top of 
that, also orbiton-orbiton interactions have to be con- 
sidered. Up to now, theoretical predictions for the two- 
orbiton contribution to <j{ui) are not available. Due to 
momentum conservation, o~(u>) is restricted to the obser- 
vation of two-orbiton processes with k to t = ki + k-2 ~ 0. 
This means that the two orbitons have opposite momenta 
hi = — k 2 with arbitrary k^, thus the orbiton dispersion 
is probed throughout the entire Brillouin zone. The to- 
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FIG. 7: (color online) Temperature dependence of — ln(T)/d 
of H0VO3 for E\\a (top panel), E\\b (middle), and E\\c 
(bottom). The band of sharp lines at 0.6 eV originates from 
crystal-field excitations within the Ho 4/ shell. 
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FIG. 8: (color online) Anisotropy of the absorption spectra 
for E || a and E \\ b for H0VO3. 



tal line width is a convolution of twice the orbiton band 
width and the width arising from vibronic coupling to 
the lattice. 

In order to test our interpretation of the feature in 
YVO3 at 0.4eV for E \\ c, we study H0VO3, which 
shows a very similar crystal structure and very similar 
magnetic and orbital ordering patterns. The absorption 
coefficient a(u) oc — In (T(u>))/d of H0VO3 shows qual- 
itatively the same mid-infrared features as observed in 
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FIG. 9: (color online) Peak frequency observed for E || a in 
YVO3 and H0VO3. 



YVO3, i.e., a peak at 0.4 eV for E \\ c in the intermedi- 
ate phase, and a peak at about 0.55 eV for E \\ a in the 
whole temperature range with the same anisotropy be- 
tween a a (oj) and <r b (oj) (see Fig. [7]). In H0VO3, the phase 
transitions occur at T s ~ 40 K and T 00 ~ 188 KP Addi- 
tionally, an absorption band consisting of several sharp 
lines is observed around 0.64 eV, which can be attributed 
to the /-/ transitions of the Ho 3+ ion ( 5 / 8 — > 5 / 7 ). The 
fine structure of this band is due to transitions between 
different Stark components. Note that the polarization 
of the incident light mainly affects the intensity of the ob- 
served lines. The close similarity of the spectra of YVO3 
and H0VO3 regarding the orbital excitations within the 
3c? shell and in particular the sensitivity to the phase 
transitions clearly show that the considered features re- 
flect intrinsic properties of the vanadates and corroborate 
our interpretation. 

What are the implications of our results for th e claim ed 
observations of orbitons by Raman scatterin g* 17 ! 18 ! 19 ! at 
energies of the order of 40 - 80 meV (see the introduc- 
tion)? The excitation energy of 0.4 eV for E \\ c implies 
that A Q + Ap w 0.4 eV. This may be consistent with a 
comparably small one-orbiton energy if the CF splitting 
on V(l) and V(2) differ substantially, e.g. A Q = 50 - 
100 meV and A^ = 300 - 350 meV. Note that a rather 
large value of A^ — A a is in agreement with the ob- 
served value of the spectral weight (see above). However, 
the spectra for E || b indicate A T < 300 meV, implying 
A r > 100 meV. 

Finally, we turn to the feature at 0.55 eV for E || a. 
Its rather high energy suggests that excitations from 
the low-lying xy orbital are involved. At all tempera- 
tures studied here, there is no significant contribution 
to ab(uj) around 0.5 eV (see Figs, [I] and [7]). This pro- 
nounced anisotropy puts severe constraints on the inter- 
pretation. Assuming an ideal perovskite structure, one 
expects <T a (ui) = (Tfc(w). In YVO3 and H0VO3, the V-O- 

V bonds are rotated within the ab plane by about 45° 
with respect to the orthorhombic a and b axes, with 
antifcrro-orbital ordering in the entire temperature range 
(see Fig. [IJ. Therefore, one expects that an exchange 
process between neighboring V sites contributes roughly 
equally to a a {uj) and <7&(u;). In contrast, next-nearest 

V neighbors are displaced parallel to the a or 6 axes. 
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One may speculate that the structural distortions, i.e., 
rotations and tilts of the octahedra, give rise to the ob- 
served anisotropy. Such a detailed theoretical analysis of 
the exchange between next-nearest neighbors including 
structural details is beyond the scope of the present pa- 
per. In order to highlight the anisotropy and to obtain 
a better view of the line shape, we plot the difference 
spectra between a and b axis of H0VO3 in Fig. [8] The 
temperature dependence of the peak frequency is shown 
in Fig. [9| Upon cooling down from 300 K, the peak po- 
sition, spectral weight and line shape all change below 
about 100 K. This is more evident in H0VO3, because 
the change is cut off at the first-order transition at T$, 
which is higher in YVO3. Possibly, this may be related 
to the Neel temperature, T N = 116K in YV0 3 and 114K 
in H0VO3P2I This points towards the entanglement of 
spin and orbital degrees of freedom, which is expected if 
exchange interactions are dominant!^ 



IV. CONCLUSION 

In conclusion, the orbital excitations in YVO3 and 
H0VO3 have been studied by optical spectroscopy. We 
focused on an absorption band observed at 0.4 eV for 
E || c. We have shown that this feature is located far 
below the Mott-Hubbard gap and that it can neither be 
interpreted in terms of phonons, magnons, weakly bound 
(Mott-Hubbard) excitons, or polaronic carriers trapped 
at impurity sites. Therefore, we identify this feature as 
an orbital excitation. However, based on the comparison 
with the data of VOC1 and with recent calculations) 15 * 16 ^ 
we have shown that this absorption peak cannot be ex- 
plained in a local crystal-field scenario, i.e., within single- 
site physics. Alternatively, we propose that this peak 
reflects collective orbital excitations, i.e., orbital excita- 
tions that are based on the exchange coupling between 
neighboring V sites. We demonstrate that the exchange 
of two orbitals between adjacent sites along the c axis in 
the intermediate phase directly contributes to a(oj). The 
energy, polarization and temperature dependence as well 
as the spectral weight of the absorption band at 0.4 eV 
are in excellent agreement with the expectations for a 
two-orbiton excitation. 

Our results call for a number of further investigations. 
Our claim can be tested directly by the direct observation 
of the dispersion with a momentum-resolving technique 
such as resonant inelastic x-ray scattering or electron en- 
ergy loss spectroscopy. Moreover, we call for theoreti- 
cal studies of the orbital exchange that realistically take 
into account the coupling to the lattice. In particular, 
a comparison of our data with predictions for the line 
shape of two-orbiton absorption is expected to reveal im- 
portant information on orbital-orbital interactions. Fi- 
nally, more detailed investigations of the exchange be- 
tween next-nearest neighbors within the ab plane are 
necessary to clarify the nature of the absorption band 
observed at 0.55 eV for E II a. 



V. APPENDIX 

Here, we discuss the derivation of the effective Hamil- 
tonian 7J or b (see Eq. |3| . Conceptually, this is even in sec- 
ond order less trivial than one might think at first glance. 
This is so because the states without double occupancy 
(states (a) and (b) in Fig. [6]) are not degenerate due to 
the differences in the energies e T and A T . This leads to 
the remarkable phenomenon that different second order 
calculations lead to different results. This stems from 
the different ways to perform the unitary transformation 
which eliminates the terms which change the number of 
double occupancies. Similar observations were made pre- 
viously in the derivation of the electron-electron attrac- 
tion mediated by phonons^. We illustrate this issue here 
by two calculations. 

Both calculations require to split the hopping part of 
Hp (second term in Eq. (lal) in two parts 



*(<£ic, 



/30 "T c q0 c /31 



h.c.) = H F+ + H F „ (10) 



with H F - = (H F+ Y and H F+ = H F++ + H F -+. The 
first plus (minus) sign indicates that a double occu- 
pancy is created (annihilated). The second sign indi- 
cates whether an electron is raised (+) or lowered (— ), 
i.e., hops from to 1 (+) of vice versa (— ). This implies 

Hp |_ = (Hp-i )t and Hp = (Hp ++ y. In detail, we 

have 

H F++ = t{c ] al c pQ h a o + 4 x c a0 n^o) (Ha) 
h f+- = t{ci c pi n al +c\ Q c al h 01 ) . (lib) 



A. Standard Unitary Transformation 

The standard approach is to determine an antiher- 
mitean operator 77 = 77+ — r\- such that 



H orh = exp(?7)i7 F exp(-?7) 



(12) 



holds. To eliminate the hopping in leading order we re- 
quire [77, Hpo] = —Hp + — Hp_ which leads to 

V++ = 7TT7^ a^ + 77, , a , a- ( 13a ) 



V+- 



U' + A a - Ae U' + A f3 + Ae 
U' - A p - Ae U' - A a + Ae 



, (13b) 



where we used Ae = ep — e a , and ?/+ = r/ ++ + r/^ as 

for the parts of the Hamiltonian, with 77_ = r/ + 77 |_, 

7] = {rj ++ y , and 77 |_ = (77^ y. In second order in 

t/0(U') we obtain ff orb = ±[77, H F+ + Hp-]. Using the 

shorthand i? rb+ = Jb^bL for the creation of orbitons we 
have to compute 

H oth+ = (1/2) ({ V++ ,H F - + } + [H F++ ,r,- + ]) . (14) 
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Explicit commutation leads to the standard result 



which implies the differential equations 



Jstan 



t 2 U' 



t 2 U' 



(U') 2 - (A Q - Ae) 2 {U'Y - (Ap + Ae) 



(15) 

For A T = = Ae this is identical to the result known 
from the derivation of the Heisenberg spin exchange as in 
(JHeisen/2)(S'+ S7 + S~ St) which implies 2 J = J Hc i se n = 



4t 2 /U'. Note that Eq. (15 1 for the exchange J becomes 
singular as soon as U' — ~>~~\A T ± Ae|. We will see that a 
smoother unitary transformation provides a less singular 
result. 



B. Continuous Unitary Transformation (CUT) 



It might surprise that the result ( 15 I is not unique. But 



we emphasize that only the matrix elements on-shell, i.e., 
without energy change, are defined independently from 
the chosen basis. All other matrix elements do depend on 
the chosen basis. Generally, a continuous change of basis 
is smoother and less singular than the one-step transfor- 
mation, see also Refs. 1561571 

The continuous change of the Hamiltonian is parame- 
terized by £ S [0, oo) and H F (£) is given by the differen- 
tial equation 



d e H F (£) = [ V (£),H F (£)}. 



(16) 



It is understood that H F (£ = 0) is given by the Hamilto- 
nian Hp in ([!]) while H F (£ = oo) is given by -fforb in 
The transformation ( 16 I shall eliminate the terms in Hp 
which change the number of double occupancies, i.e., the 
kinetic part Hp + + Hp- . Hence we parameterize 

H F++ (£) = A 1 (e)cl 1 c go n aQ + B^c^c^hgo (17a) 
H F+ -(£) = A (e)cl a c gi h al +B (i)cl Q c al h 31 , (17b) 

while Hpq remains constant in linear order in t; the 
operators H F —{1) = (H F++ (£))^ and Hp_ + {£) = 

(Hp^ (£)y follow by hermitean conjugation. 

The crucial choice is the one for the infinitesimal gen- 
erator r){l). Our aim is to eliminate processes which cre- 
ate or annihilate excitations of the order of U' . Such 
an elimination can most easily b e done by the Miclke- 
Knetter-Uhrig generator t/mku (l^ 58 ! 59 ! 60 ! which consists 
of the terms in the Hamiltonian increasing the number 
of excitations and of the negative terms in the Hamilto- 
nian decreasing the number of excitations 



VmkvW :=H F+ (l)-H F _(£), 



(18) 



for a general discussion see also Ref . [52 With this choice 
one obtains in linear order in t 



d e Ai = -(U' + A a - As)A x (20a) 

d e A = —(U' — Ap— As)A (20b) 

d t B x = -(U' + A l3 + Ae)B 1 (20c) 

d e B = -(W - A a + Ae)B . (20d) 

The solutions consist in decreasing exponential functions 
starting at t for I = because we assume all the energy 
differences in the parentheses in (20 1 to be positive, i.e., 
U' dominates the other energies, see Eq. [3] 

The orbital exchange is obtained by equating the sec- 
ond order terms in t in (16) which implies 



d t J{£)b% = 2{H F++ (e),H F _ + (£)}. (21) 



Since the right hand side is given by the solutions of ( 20 1 
an integration suffices to provide Jcut = J{£ 



oo 



Jcut = 2 / {A^i) A (£) + B^i) B (£))de (22a.) 
Jo 

_ t 2 t 2 

U' -S + ecE + U' + 6- ecE 

(22c) 



(22b) 



(U') 2 -(5- ecEf 



where we use the shorthand 5 for the crystal-field levels 
(see Eq.[6|, and E denotes the applied electric field. This 
is the result used in the main part of the article, see Eq. 

13 

Note that we retrieve the well-known result for the 
Heisenberg exchange of 2 J = JHcisen = 4i 2 /J7' for E = 
and (5 = 0, i.e., equivalent V sites. 

Even more interesting is that Jqut 7^ Jstan- In partic- 
ular, the individual excitation energies A T do not occur 
in Jcut in ( 22c ) but only their difference (see Eq. |6| . 
Hence a regime exists with 5 = and A T — > U' where 
Jstan diverges while Jcut remains unaffected. So the 
CUT result is less singular. Moreover, the expression for 
Jcut is simpler than the one for J s t a n- 

Tracing back from where the difference between the 
standard and the CUT results originates we have to com- 
pare Eqs. ( 14) and (21 1. In the standard calculation ( 14 1 



there is a striking asymmetry between the two operators 
which are commutated. Only one of them (77) carries in- 
formation on the excitation energies. In the CUT calcu- 
lation (21 1 both commutated operators carry the depen- 



dence on the excitation energies in the same way by their 
dependence on £. This implies also that a single commu- 
tation suffices because two commutators are equal while 



there are two different ones in (14) 



d e H F+ (£) = -[H F+ (£),H F0 } 



(19) 



For all the above reasons we favor the CUT derivation. 
We stress, however, that in the regime relevant for YVO3 
the difference between J s t a n and Jcut is quantitatively 
of minor importance. 
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